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A sequence of bursts observed in an intermittent time series may be caused by a single avalanche, 
even though these bursts appear as distinct events when noise and/or instrument resolution impose 
a detection threshold. In the Bak-Tang-Wiesenfeld sandpile, the statistics of quiet times between 
bursts switches from Poissonian to scale invariant on raising the threshold for detecting instantaneous 
activity, since each zero-threshold avalanche breaks into a hierarchy of correlated bursts. Calibrating 
the model with the time resolution of GOES data, qualitative agreement with the inter-occurrence 
time statistics of solar flares at different intensity thresholds is found. 
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A simple picture of intermittency in turbulent flows 
places rare, hot regions that dissipate energy inside a 
cold laminar sea. As the hot regions evolve, they main- 
tain a clustered structure, and the dissipation, occurring 
at small scales, remains correlated at large ones. For 
instance, examining the sky at night, one sees stars, 
galaxies and clusters of galaxies against a dark back- 
ground 0,0. Many striking examples of turbulent inter- 
mittency occur in astrophysical, space or laboratory plas- 
mas, such as flaring events in the solar corona 00, mag- 
netic substorms 0, bursty bulk flows M, auroral emis- 
sions , turbulence in the solar wind [1, , or bursts 
observed in RFX E3 and RFP experiments [l2| . 
Important parameters include the Reynolds number (s), 
R = VIq/v (or R m = Vlo/rj, etc.), where v is the viscos- 
ity, r\ is the magnetic diffusivity, and V is the velocity dif- 
ference over the integral scale Iq. For increasing R, char- 
acteristic widths of the dissipating regions decrease, while 
their intensity increases. Thus turbulence becomes an 
"on/off" phenomena as the appropriate Reynolds num- 
ber's) become large 0EH- 

The controversial hypothesis that turbulent intermit- 
tency may be a manifestation of self-organized criticality 
(S OC) has been discussed by Bak and others 0, 0, IE 
Il4 Il5l ll6L Il7l Il8l Il9l 1201 |2l| . In this scenario, intermit- 
tent energy dissipation is a stick-slip or threshold process. 
Each slip can trigger further slips - either through short 
or long range interactions. Eventually, a regime materi- 
alizes where sparse, sporadic avalanches of intense energy 
dissipation interrupt laminar regions of space-time which 
rest near static equilibrium with little dissipation despite 
the continuous, global input of energy. 

Up to now, SOC has mostly been studied with vividly 
plain models, such as the original Bak-Tang-Wiesenfeld 
sandpile (BTW), that ignore many features of hydrody- 
namic, plasma or other kinds of turbulence. Rather, they 
generate patterns of rapid energy dissipation in space 
and time, which is the hallmark of intermittency. On 



the other hand, most studies of turbulence [13| examine 
structure functions and multiscaling phenomena [22! I"2.'"{l| , 
and/or nonlinear instabilities and coherent structures, 
etc. 0, all of which are associated with inter- 
mittency but do not directly characterize the bursts of 
energy dissipation in space and time. As a result, some 
comparisons that have been made are superficial and 
should be made more definitive. Further, one argument 
used so far to distinguish SOC from turbulence is mis- 
leading and erroneous. 

While BTW and other SOC models exhibit a broad 
distribution of avalanche sizes and durations, which are 
comparable to e.g. solar flare data, a marked differ- 
ence has been noted regarding the time intervals be- 
tween bursts. For instance, Boffetta et al. 26] found 
that the distribution of times between flares exhibits 
power law statistics, while intervals between subsequent 
avalanches in BTW are (approxim ately) P oissonian. Fur- 
ther, they and other groups E3, El 111 E3 found that 
shell or reduced magnetohydrodynamics (MHD) models 
gave a better description of the waiting time statistics, 
and some E3 EH HE E3 used this to distinguish SOC 
from turbulence, or to question the applicability of the 
SOC paradigm for magnetically confined plasmas in ther- 
monuclear research [lfj . 

We make three key points. First, intermittent bursts 
can never be detected, nor distinguished from the back- 
ground, at arbitrarily low thresholds. For instance, the 
studies above comparing reduced MHD or shell models 
with flares or bursts in man-made plasmas use a thresh- 
old for defining bursts. Such a threshold is realistic be- 
cause the emission associated with e.g. flares decays 
slowly after a local peak, allowing overlaps with subse- 
quent peaks. Although a threshold is unavoidably con- 
nected to the precise definition of events, robust features 
may be observed with rescaled distributions measured at 
different thresholds [H|2^|. For such time series, event 
durations and quiet times between them are measured on 
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FIG. 1: (Color online) Time series of the instantaneous num- 
ber of topplings, n(t), in the two dimensional BTW with 
L = 2048 and AT = 100. Dark (red online) areas and the 
black line represent the actual signal and the signal with some 
additive noise, respectively. Shaded regions indicate a detec- 
tion thresholds n c = 10 (yellow, dashed line) and n c — 100 
(blue, dot-dashed line); any signal below n c is considered un- 
detectable. Intervals between a rise of the signal above n c 
and its next fall below constitute measured durations td of 
events, followed by a "quiet time" t q until the next rise. The 
"waiting time" between two consecutive rises is t w = t d + t q . 

a single clock with equal precision. Thus one can consider 
the hypothesis that the sequence of bursts arises from 
a single avalanche observed at finite detection threshold. 
The bursts within an avalanche in a SOC system can be 
expected to be correlated in space and time, being part 
of the same, critical process. 

To test this idea, we study a BTW sandpile |b| in 
which the time unit is that of a parallel update step, 
and consider both slow driving (A) and running sand- 
pile j^jj, H3, E01 (B) conditions. In the time series of the 
activity, n(t), bursts are defined as consecutive intervals 
during which n(t) > n c , where the detection threshold for 
events is n c > (see Fig. A data analysis technique 
analogous to that recently developed by Baiesi et al. |28j 
to examine inter-occurrence statistics of flares provides a 
direct comparison between results from numerical simu- 
lations of BTW and solar flares. It turns out that BTW 
shares several features with the intermittent statistics of 
flares, as shown below. 

The BTW sandpile consists of a L x L lattice with a 
discrete number Zi of sand grains occupying each site i. 
We study two versions: (A) in the slow driving limit, a 
grain is added to the pile at a randomly chosen location 
when the previous (n c = 0) avalanche ends. The du- 
rations (td) and quiet times (t q ) then refer to intervals 
between local peaks within each avalanche and statistics 
are obtained over many avalanches. In the second case 
(B) , one grain of sand is dropped every AT update steps 
at a randomly chosen site. In both cases, at each up- 
date step, t, all sites that exceed a threshold for stability, 



Zi > z c — 3, topple in parallel by distributing a single 
grain of sand to each of their four nearest neighbors or, 
for boundary sites, over the edge of the lattice. Taken 
as the instantaneous dissipation signal, the activity n(t) 
is the number of unstable sites toppling at each parallel 
update step. In model B, if AT > (tdjL, the sequence 
of topplings is also interrupted by instances where the 
activity completely stops (n(t) = 0) 31]. The quantity 
(td)L is the average duration of avalanches on a lattice 
of scale L in the stationary state of model A. For both 
A and B, consecutive stopping points separated by in- 
tervals where n(t) > delimit n c = avalanches. The 
time series of Model B has similar character to the solar 
flare data studied in Ref. [2^], with a broad distribution 
of events that exceed each threshold n c , albeit with a 
finite-size cutoff curtailing the power-law tail observed in 
Fig. 2a of Ref. 0. 

Consider an observer who measures the global activ- 
ity sequence with a finite error, so that the time series 
she records is n b s (i) = n(t) + rj(t). For instance, let 
r](t) be an independent random number uniformly dis- 
tributed between and 15. The effect of this noise is 
shown in Fig. ^for model B. One way for the observer 
to separate the signal from the noise is to increase her 
threshold for detecting events, and coarse-grain her unit 
of measurement. In observing natural phenomena, such 
as flares, these detection thresholds are an intrinsic and 
unavoidable part of the measurement. Hence, we study 
the original time series together with a finite threshold 
n c > (e. g. n c — 100 as in Fig.^l to distinguish bursts, 
considering all instances with n(t) < n c to have no activ- 
ity. For sufficiently large n c , the event statistics that the 
observer measures (at large times) are the same as the 
actual statistics without noise. 

As Fig. [21 shows, on increasing the threshold n c from 
zero, the distribution of quiet times t q for model B 
switches from an (approximately) exponential distribu- 
tion to a power law, 

P qu iet(t g )~t-^ with 7™ = 1.67 ±0.05. (1) 

For model A one observes an even cleaner and broader 
scaling regime, with the same r ) q 3TW |4l|. In model A, 
the power law tail arises from the correlations of bursts 
within each avalanche. Further, the Abelian property of 
BTW assures that the power law behavior of quiet times 
in model B cannot be due to overlapping avalanches. 
Thus, the natural introduction of detection thresholds 
leads to the discovery of a hierarchical sequence of corre- 
lated bursts (or sub-avalanches) within a large avalanche. 

A similar data analysis by Baiesi et al. [28j for flares, 
detected as intervals during which the emission intensity 
in the GOES time series exceeds /, used a particularly 
simple scaling ansatz for the quiet time distribution, 

Pquiet(tq\I) = J-^-f flare ( jA-) ■ ( 2 ) 
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FIG. 2: (Color online) Distribution of quiet times t q in BTW 
for different thresholds n c for model A and model B with 
AT = 100 (groups of curves are offset vertically). In model 
B the distribution changes from approximately exponential 
at n c = to a power law for increasing n c . The straight 
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reference lines decay as t d 1 
Eq. p. 



suggesting 7 f J w = 1.67(5) in 
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FIG. 3: (Color online) Rescaled distribution of quiet times in 
BTW with L = 2048, AT = 100 for different thresholds n c 
(black lines, open symbols) and for the solar flare data cor- 
responding to "ALL" at different thresholds I in Fig. 3(b) of 
Ref. |2gj| (colored lines, full symbols). The quiet times t q have 
been divided by the average quiet time (t q ) at each thresh- 
old, and the distributions rescaled to preserve normalization. 
BTW and flare data are similar in the intermediate regime. 
The former are shifted up by one unit on the log scale. 



Here ffiare is a scaling function and the average quiet 
time for a given threshold, (t q )i, provides a factor that 
collapses distributions measured at different /. While 
BTW does not obey this particular scaling ansatz, for 
comparison we apply Eq. J2J) to BTW with different val- 
ues of n c . This enables us to study the dependence of 
two dimensionless quantities {t q )P qu iet vs. t q /{t q ) at dif- 
ferent thresholds for both flares and BTW. As shown in 
Fig. we find that although the scaling functions f(x) 
are clearly different for x <C 1, they are similar for inter- 
mediate arguments. 

Event durations, td, observed at different thresholds in 
BTW and for flares occurring during solar minimum, are 
shown in Fig. 0] 42] . In order to compare with BTW, 
we set the time resolution of the GOES data equal to 
the time resolution of BTW. One minute is set equal to 
one sweep. Using this straightforward calibration, the 
statistics of event durations are similar as well. In both 
cases, power law behavior is observed, Pdur{td) ~ t^ ldur . 
For flares at solar minimum jdur = 2.0 ± 0.1 [2^, while 
this exponent the n c 1 bursts in BTW is smaller. 

The inset of Fig. 0] shows how the apparent power law 
behavior for burst durations in BTW changes with in- 
creasing threshold n c . At sufficiently large n c a plateau 
appears in the function t d Pdur{td)- This suggests a crit- 
ical exponent at high thresholds, Jdm- W = 1-67 ± 0.05, 
similar to j q TW ■ The figure also demonstrates that the 
scaling behavior of burst durations, for the system sizes 
that were studied, differs at large n c from that at n c = 0. 
In the latter case, Liibeck and Usadel [13] determined a 
critical exponent 7^ w/ (" c = 0) ~ 1.48, while Stella and 
De Menech [33| found multiscaling. 
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FIG. 4: (Color online) The distribution of event durations 
for BTW model B with L = 2048, AT = 100 and for flare 
data corresponding to "minimum" at different thresholds I 
in Fig. 4 of Ref. _28j . One parallel update step in BTW has 
been calibrated to the time resolution of the GOES data, or 
one minute. Although the critical exponents are different, 
the overall distributions are roughly comparable. The inset 
shows that, for increasing n c , the behavior changes from that 
in Refs. |32l |33| towards a plateau indicating an asymptotic 
exponent m 5/3 for BTW at large thresholds [n c : (•), 15 
(+), 255 (o), and 1023 (V)]. Similar behavior is found for 
model A (not shown). 

If the sequence of quiet times (t q )i are uncorrelated, 
the cumulative variable yi(j) — Yll=j (tq)i exhibits 
diffusive behavior, with an average variance scaling as 
g = ((y 2 ) - (y) 2 ) 1/2 - l H , with H = 1/2 34]. Our mea- 
surements of this quantity confirm that quiet time inter- 
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vals for BTW at different thresholds are uncorrelated, in 
all cases giving H = 1/2. However, at solar minimum, 
we get H ~ 0.62 (see also Therefore, BTW does 

not reproduce correlations between quiet times for flares. 

In fact, a variety of SOC models exhibit power laws 
in times between events. Previous analyses |3(| IU H|| 
using a detection threshold considered a completely dif- 
ferent limit than that discussed here, namely that of an 
infinite time scale separation between the driving rate 
and the durations of events. Consequently, a threshold 
was imposed on a variable related to the total size or du- 
ration of each n c — avalanche, rather than its instan- 
taneous dissipation. However, as explained earlier our 
limit of overlapping time scales for durations and quiet 
times may be the correct one to describe intermittency in 
turbulence. Indeed, cellular automata models of labora- 
tory plasmas are running sandpiles; see e.g. Newman et 
al |33|. Within this scheme, some works |2oLl4*c| . also us- 
ing a threshold, have found power law quiet times for the 
Hwa-Kardar |3f| running sandpile driven at a sufficiently 
high rate. This was claimed to be due to interactions be- 
tween overlapping avalanches. In contrast, here we show 
that the BTW model, considering one avalanche at a 
time, generates a power law distribution of quiet times 
when a finite detection threshold is used. Its Abelian 
property assures that this correlation remains the same 
when the model is driven at a finite rate, in agreement 
with the results shown here. It is unlikely, though, that 
the Abelian property of BTW is essential to getting a 
scale free distribution of quiet times, although this re- 
mains to be clarified by studying other models in a sim- 
ilar way. 

In conclusion, we have demonstrated that SOC re- 
mains a viable alternative for the explanation of intermit- 
tent dissipation in turbulence by comparing results from 
numerical simulations of a BTW sandpile with solar flare 
statistics. By including an inevitable detection threshold 
in the analysis of BTW, as well as allowing an overlap 
of time scales between burst durations and quiet times, 
qualitative correspondence is obtained for the power law 
statistics of inter-occurrence times for solar flares. Stud- 
ies of more physically realistic SOC including detection 
thresholds for short time dissipation in the whole system 
may improve quantitative agreement. 
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